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ABSTRACT 



This thesis proposes and investigates a new approach in attacking the 
"need to calibrate" problem of the shipboard HFDF systems. It is based on 
measuring the system response to multiple onboard near-field sources. The test 
is performed along with the standard calibration, and the antenna responses to 
the near-field sources are stored in a reference database. Whenever a 
modification is made to the topside, the near-field test is repeated and the new 
results are compared to the near-field reference database. A significant 
difference may indicate a need to perform a full system calibration. The 
calibration procedure was simulated using the numerical electromagnetics code 
PATCH. Calculations show that the antenna responses for near and far-field 
sources are comparable when topside changes were introduced. 
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I. INTRODUCTION 



"Thalassocracy" is a word introduced to Western civilization by its earliest 
historians whose works survive, Herodotus and Thoucydides. "Thalassa," in 
Greek, means the "sea." To these men and their contemporaries, thalassocracy 
was a concept which meant, loosely, "maritime supremacy," i.e., the control of 
the sea lanes and islands by one state to insure its territorial integrity and thus 
its economical prosperity. "Thalassocracy," therefore, is a universal concept, 
with technology being the corner stone. Technology has been a driving force in 
defense since the man threw the first stone thousands of centuries ago. As 
weapon systems advance in complexity, effectiveness, and in the ability to 
evade countermeasures, the user demands more of the early warning detection 
and location systems that support their efforts. 

Electronic warfare (EW) is a distinct and well-defined major function in na- 
val operations. In its broadest sense, EW encompasses the employment of all 
devices and equipment that radiate or receive electromagnetic energy. The in- 
herent weakness of an electromagnetic signal is its susceptibility to detection 
and interference because of its free-space propagation characteristics. It is the 
exploitation of this weakness which accounts for a major expenditure of effort in 
the electronic warfare arena. 
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Of particular interest in this thesis is the field of tactical radio direction find- 
ing for early warning and position location of targets. High frequency direction 
finding (HFDF) has been widely used since World War II. Tactical commanders 
at sea had their first look at the utility of HFDF during the battle of the Atlantic. 
HFDF stations were able to provide relatively accurate locations of German U- 
boats, enabling seaborne units to successfully intercept the submarines. HFDF 
was eventually installed on board naval surface ships. With the advent of 
cruise missiles and the need for over the horizon targeting (OTHT), shipboard 
HFDF was seen as a means of obtaining targeting information without divulging 
one’s own position. Passive signal intelligence offers the tactical commander a 
viewpoint into the intentions of his adversaries; HFDF further enhances battle- 
field effectiveness by providing knowledge of the enemy’s location. 

One of the basic steps in the development of an HFDF system, or any other 
system, is the determination of the calibration method needed to ensure its op- 
eration within the limits set by either the manufacturer or the user. The scope of 
this thesis is to investigate a method of evaluating the response of the antenna 
elements of such a system, both before and after modifications to the ship, thus 
providing an indication of the need to calibrate. This method is based on the ex- 
citation of the system to multiple near-field sources, and will subsequently be re- 
ferred to as the "near-field test method." 
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Chapter II is an overview of HFDF systems and their capabilities. In addi- 
tion, the Classic Outboard and Combat DF systems, the associated hardware, 
and their operation are briefly described. 

In Chapter III the present method of calibrating antenna systems and its 
applicability to the HFDF systems are discussed. This method requires that the 
ship circle a transmitter while the system response is measured. Currently there 
is no reliable indicator as to whether a topside modification will affect the system 
DF accuracy. Therefore, the system is recalibrated after any significant modifica- 
tion. This approach is effective but very costly and time consuming. 

In Chapter IV the proposed near-field measurement approach is described 
and simulated using the electromagnetic code PATCH. A numerical model of 
the DD 963 (Spruance Class Destroyer) is used to compute the deck-edge an- 
tenna response to a range of near-field sources. It is shown that the responses 
are sensitive to topside modifications if a complete set of transmit locations is 
considered. 

Finally, in Chapter V, the data is summarized and conclusions and recom- 
mendations for future work are presented. 
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II. OVERVIEW OF HFDF METHODS 



A. HIGH FREQUENCY DIRECTION FINDING (HFDF) 

Determining the location of a remote transmitter is important to the Navy in 
search and rescue, in forecasting the actions and intentions of a potential en- 
emy, and in targeting "fire-and-forget" missiles. There are several methods for 
deriving geolocation information at HF, and various algorithms attempt to exploit 
attributes of the schemes. Multiple sites can be used to take advantage of the 
signal time-difference-of-arrival (TDOA) or bearing angle-of-arrival (AOA). Per- 
haps the most versatile class of methods are those in the category of the single- 
station-location (SSL) such as the Classic Outboard and Combat DF systems. 
SSL refers to methods by which a signal source is located from information at a 
single receiving site such as a ship. Clearly this requires that some information 
or assumptions about the propagation path through the ionosphere (sky wave) 
and over the earth's surface (ground wave). 

A HFDF system serves several purposes: 

1. integrate shipboard tactical location intelligence with longer range, 
shore-generated bearing information, 

2. exploit the strength of the system while avoiding its weaknesses, and 

3. provide this location information in a timely manner to users. 
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Position location of high frequency radio transmitters consists of three se- 
quential operations: 

1. detection of the signal, 

2. determination of the direction of emitter (DOE), and 

3. computation of the position estimate. 

Errors and inaccuracies can be introduced during any of these operations, 
and if left uncorrected will result in position estimates with little tactical value. 
The remainder of this section will briefly explain how errors are introduced in 
each phase of position locating. 

The first phase, detection of an HF radio transmitter, is usually thought of 
as a binary operation: it either occurs or does not occur. However, under certain 
signal-to-noise (S/N) conditions detection may seem to have occurred when in 
fact it hasn't, or it may not have occurred when in fact it has. This contradiction 
is caused primarily by two factors: noise and high signal density (interference). 
The lack of an observed detection when in fact the radio signal is present is nor- 
mally caused by noise at the HFDF site. The noise can be self-generated at the 
site or come from outside sources that cause interference. 

The detection of a radio signal at the frequency of interest does not guaran- 
tee that the transmitter of interest was detected. The HF spectrum is crowded 
with multiple users employing the similar signals on the same or near to the 
same frequency. Multiple HFDF sites may report detections on a frequency and 
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in fact not be hearing the same transmitter. Detection errors are most likely to 
occur when there does not exist an HF propagation path between the transmitter 
of interest and one or more HFDF sites. 

The second phase of position location is determination of the direction of 
emitter (DOE). DOE is defined as the azimuth angle from the HFDF site in the 
direction of the transmitter of interest. After detection, the HFDF site reports the 
direction-of-arrival (DOA) of the transmitted signal. When the reported DOA dif- 
fers from the true DOE an error has been introduced. The true direction of the 
transmitter is never known precisely in practice, and thus only an estimate ex- 
ists ( which will be referred to as bearing). 

Large bearing errors are usually the result of detection errors; smaller more 
insidious bearing errors are the result of a variety of factors. The most prevalent 
is the error caused by propagation irregularities. Long range HF communication 
is made possible by the ionospheric layers which act as reflecting mediums 
propagating HF signals around the curvature of the earth as shown in Figure 1. 
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The ionospheric layers are not true reflectors, but are instead refractors 
that bend the signal back in the direction of the earth's surface as shown in Fig- 
ure 2. As the signal is refracted the wavefront can be distorted away from its 
direction of propagation prior to refraction. (The plane defined by the transmitter, 
receiver, and reflection points illustrated in Figure 2 may not contain the earth's 
center.) The magnitude of the distortion will determine the degree of the bearing 
error. 



• Apparent reflection point 




Figure 2. Virtual height of an ionospheric reflection . 



Bearings can also be different than the true DOE if reflecting and scattering 
objects exist in the vicinity of the HFDF antenna. Site errors can be reduced with 
careful site selection and empirically developed corrections (background sub- 
straction methods). 
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The final stage of the position locating problem is the actual computation of 
a position estimate or "fix." In the absence of errors, the position estimate and 
the transmitter of interest are collocated and for a single station can be deter- 
mined by phase and time-of-arrival differences between adjacent antennas. 

Fix algorithms have been developed to solve the estimation problem with a 
variety of mathematical techniques. Ideally, the fix algorithm must be capable of 

1. differentiating good bearings from bad, 

2. adjusting good bearings to correct for small bearing errors, and 

3. providing an indication of the estimates accuracy. 



B. OVERVIEW OF CLASSIC OUTBOARD AND COMBAT DF SYSTEMS 

The Classic Outboard and the Combat DF systems are installed on most 
Spruance class destroyers (DD 963) as an integral part of their EW capability. 

A block diagram of Classic Outboard is shown in Figure 3. 

Both of these systems employ an array of shielded half-loop antennas with 
ferrite cores (Figure 4). There are a total of 24 antennas mounted on the surface 
of the ship and spaced around the perimeter of the deck as shown in Figure 5. 
They operate from 4 to 30 MHz and are referred to as "deck-edge antennas." 
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Figure 3. Block diagram of the Classic Outboard System. 
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Figure 4. Ferrite loaded loop antenna. 
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Figure 5. Antennas locations on the ship. 
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1. Deck-edge Antennas 

These deck-edge antennas are attached horizontally (in most cases) 
directly to the skin of the ship as shown in Figure 6. They serve as current 
probes that sense the vertical component of the current induced on the ship's 
surface. This current is primarily due to a vertically polarized incident plane 
wave, but can also be induced by horizontally polarized waves due to ship 
scattering. 



Antenna 

shield 




Figure 6. Ferrite loop on the ship's surface. 

Current probes (sensors) have been used for many years in different areas 
of technology [Ref.1], They are devices which convert an electromagnetic quan- 
tity of interest into a voltage or current at some terminal points. A design proce- 
dure introduced by Pettengill, et al [Ref. 2] provides an equation which relates 
the minimum required volume and the optimum dimensions for a fully wound fer- 
rite loaded receiving antenna within the range of 0.5 to 30 MHz 
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where, 



W = antenna volume 



F r = receiver noise factor 



k = Boltzman constant 



B = bandwidth 



T = absolute temperature 



E mi „ = minimum induced emf 



p 0 = air permeability 



fip = ferrite core permeability 



f= operating frequency 



S/N = signal-to-noise ratio 



The important feature of this equation is that the signal-to-noise ratio (SNR) is 
roughly proportional to the operating frequency. 

The probe is calibrated by the manufacturer during production by simply 
passing a current of known magnitude and frequency through the probe and 
measuring the voltage at its terminals. The result is a set of calibration curves 
that relate the voltage (V) to the current (I) and the transfer impedance (the ratio 
Z=V/I). 

2. Direction Finding Algorithm 



reduce the error in the bearing measurement. One is to use widely dispersed an- 
tenna locations around the perimeter of the ship, thus providing a wider aperture 
for phase measurement. The exact placement of the antennas is determined 



Outboard and Combat DF incorporate several techniques in order to 
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empirically to minimize the interactions and interference from objects and struc- 
tures on the ship. 



Secondly, rather than using an additive or subtractive approach to finding 
the DOA, the emitter direction is isolated using a technique called correlation in- 
terferometer direction finding (CIDF) . The direction is obtained using two sets of 
data. One consists of the phase and amplitude voltages which result from a 
known transmitter at a given frequency and relative bearing. This is referred to 
as the "reference set," and its measurement process as "calibration." The se- 
cond set consists of phase and amplitude voltage measurements received by 
each of the antenna array elements from the HF signal of unknown azimuth (the 
"measurement set"). The azimuth at which maximum correlation between the two 
sets of data occurs is taken to be the DOA of the incoming signal. 

The degree by which the measured antenna voltages matches each of the 
test voltages is called the complex correlation coefficient (C) and is computed by 
the formula [Ref. 3] 



Ln=l[^4>/W(<h)] 



( 2 ) 
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where, 



U n ((. t>j) = complex measured voltage for antenna n, 

V/ n (<i) ( ) = complex stored voltage for antenna n, 

(jjj = radiation angle of incidence, 

4> t = angle of stored voltage set, 

and the asterisk denotes complex conjugate. Note that C is normalized so that 
0 <_ C_< 1. For Classic Outboard, only vertically polarized transmissions are 
stored, while the Combat DF system uses both vertically and horizontally polar- 
ized transmissions to develop the reference database. 

To obtain a DOA, the phases and amplitudes of an unknown transmission 
are measured, and the resulting set of measurements is compared with the 
stored patterns for best fit. The closest angle-of-arrival estimate is taken as the 
one with the highest C value by polynomial interpolation, as illustrated in Figure 
7 . 
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Single Angle of Arrival (SAOA) 



Figure 7. Polynomial interpolation to find SAOA. 

The appeal of the CIDF technique is that shipboard scattering effects are 
automatically included, at least in principle. If a plane wave is incident from the 
direction ^ , then it should excite the exact same set of voltages as the reference 
signal from the same direction. Thus <}>,= ^ and C=1. In practice, does 

not always result in perfect correlation because of environmental differences be- 
tween the operational measurement and the baseline measurement. In other 
words, propagation conditions are variable due to ocean and ionospheric condi- 
tions. These cannot be controlled and therefore always introduce some uncer- 
tainty into the DOA measurement. 

In addition to environmental variations, any significant change to the ship 
between calibration and measurement will also affect accuracy. Precisely which 
change is significant depends on the angle of arrival of the signal and the loca- 
tion of the change. Clearly, a minor change at the bow will not affect antenna 
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responses astern, but can certainly affect forward antenna performance. The net 
effect on DOA estimation depends on the relative importance of each antenna's 
response in computing C. For instance, if the magnitude of an antenna's com- 
plex voltage is small relative to other antennas, then large phase errors will have 
little effect on the DOA estimate. 

The "need to calibrate" procedure introduced in Chapter III is a means of 
evaluating the accuracy of the calibration data-base after a modification is per- 
formed on the topside. Thus, it provides a solid criterion to decide whether a 
need for full-scale calibration of the system exists. 
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III. CALIBRATION METHOD 



A. PRESENT METHOD OF CALIBRATION 

One of the basic steps in developing and maintaining a HFDF system is its 
calibration and performance evaluation. The present method of calibrating HFDF 
shipboard systems is to circle the ship about an anchored buoy while monitoring 
signals from a shore site as depicted in Figure 8. 




R=10nm 

Figure 8. Calibration geometry. 

The procedure requires that the ship circle offshore while a transmitter 
frequency sweeps over the band of interest. At a discrete set of frequencies the 
DF system response is recorded along with the ship's bearing. Together they 
completely provide the reference voltage set V(<j> ( ). However, because of ship 
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pitch and roll and the fact that the course is not always circular, the apparent 4> t 
is slightly different than the actual value. 

A complete calibration requires several days and as a consequence is very 
expensive. Currently the calibration is repeated after every "significant" modifica- 
tion, with significant being determined by intuition and past experience. The 
method proposed in the next section attempts to quantify the need to calibrate. 

B. THE PROPOSED METHOD 

This thesis proposes and investigates a new approach in attacking the 
"need to calibrate" problem of the HFDF system. It is based on measuring the 
system response to multiple onboard near-field sources. The test is performed 
along with the standard calibration described in the last section, and the antenna 
responses to the near-field sources are stored in a reference data base. When- 
ever a modification is made to the topside, the near-field test is repeated and the 
new results are compared to the near-field reference data base. A significant dif- 
ference may indicate a need to perform a full system calibration. These mea- 
surements may provide additional information beyond the "need to calibrate." 
For instance, they can be used to compare the effect of identical modifications 
on two different ships. If the near-field test results yield the same change for 
both ships, recalibration might need to be performed only on one ship; changes 
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to the calibration data for the second ship could be inferred by comparing it's 
near-field data with the first one's. 

In order to validate the proposed approach, a simulation was performed 
which is briefly described as a three step procedure: 

1. compute the response of the antennas due to a source placed at various 
locations on the ship within the near-field of the antennas, 

2. perform the same set of computations after having introduced a simulated 
modification to the topside of the original model, and 

3. compare the results to determine whether any difference is significant. 

There are several near-field parameters that are varied in the simulation to 
determine their effect on the outcome and, perhaps, arrive at some optimum or 
acceptable values. They include: 

1. number and location of the transmitters, 

2. transmitter power level, 

3. source proximity to the ship modification, 

4. source proximity to antennas, and 

5. type and severity of the modification. 
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IV. NUMERICAL MODELING AND SIMULATION OF THE CALIBRATION 



METHOD 

A. NUMERICAL MODELING OF THE SHIP 

The "need to calibrate" can be demonstrated in either of three ways: 

1. full scale measurement, 

2. scaled-model measurement, or 

3. computer simulation. 

The last, computer simulation, is the easiest because it does not require a ship 
(or model) and test instrumentation. Changes can be quickly simulated on the 
computer, whereas corresponding changes on a ship would be extremely costly 
On the other hand, the numerical model does not have the fidelity of the full- 
scale model. However, this is not a critical problem because the "need to cali- 
brate" procedure is only looking for relative changes in the antenna responses, 
not absolute values. 

The ship modeling as well as the procedure simulation was performed us- 
ing the code PATCH. PATCH is a frequency domain electromagnetic code 
based on a method-of-moments (MoM) solution to the electric field integral 
equation (EFIE), and can be applied to objects of arbitrary shape, both open 
and closed. 
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The code has many capabilities and they are fully described in the User's Manu- 
al [Ref. 4], However, it is convenient to highlight the ones which were exercised 
during the course of the study. 

PATCH uses planar triangular patches for modeling complex objects be- 
cause they can easily conform to surfaces of general shape and allow variable 
densities over the surface of the object. The model of the DD-963 (Spruance 
Class Destroyer) is shown in Figure 9. The lines denote edges of flat triangular 
patches. The requirement for patch dimensions on smooth surfaces is that 
the edge lengths should typically be no longer than 1/5 to 1/10 of a wavelength. 
This easily satisfied at frequencies below 10MHz where the wavelength is 30m. 




Figure 9. Patch model of the Spruance class destroyer (DD-963). 
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The code makes full use of the existing symmetries in the geometry of the 
problem by placing perfect electric or perfect magnetic conductors in place of the 
symmetry planes ( x=0, y=0, z=0). Thus, the surface of the sea is represented by 
a perfect electric conductor placed at the z=0 plane. 

The near-field source is a very thin tapelike structure which simulates a 
monopole on the ship deck. The rule of thumb for modeling thin wires is that a 
circular wire of radius r is equivalent to a tape width of 4r. The simulation of the 
near-field procedure used the following data: 

1 . frequency: f = 6 MHz, 

2. monopole driving voltage. V=100 V, 

3. monopole antenna length : X / 10, and 

4. monopole antenna width : X / 30. 

The deck-edge antennas of the HFDF system are assumed to be located 
along the edge of a patch as shown in Figure 10. The code provides the current 




Figure 10. Antenna location on a patch edge. 



25 



densities perpendicular to each edge for all patches. The selected antenna ele- 
ments are located along horizontal edges because the current crossing them is 
due to the vertical component of the incident electric field. The current density 
J, is constant along the length of an edge. Thus, if the deck-edge antenna 
length is 0.5m and the edge is L, the current sensed by the loop is 1= J s L. 

The locations considered in the simulation are: 

1. Bow : Edge number 10, B-10, Figure 11. 

2. Midships : Edge number 1762, M-1760, Figure 12. 

3. Astern : Edge number 1 06, A-1 06, Figure 1 3. 



Edge number B-10 
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Figure 11. Antenna element on edge number 10 (bow). 



26 




Figure 12. Antenna element on edge number 1760 (midships). 



Edge number A-1 06 
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Figure 13. Antenna element on edge number 106 (astern). 
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B. NEAR-FIELD MEASUREMENT SIMULATION 



The first step is to excite the antenna elements with a source located at 
various positions on the surface of the ship. PATCH is used to compute the cur- 
rent densities on the ship due to the source. The antenna responses are deter- 
mined from the current passing across the antenna length. The coordinates of 
the sources considered in the simulation are shown in Figure 14. The coordi- 
nates of the strip are as shown in Table 1 : 



TABLE 1. STRIP COORDINATES. 



LOCATION 


*1 ■ *4 


* 2' X 3 


y 


Z 1 - Z 4 


Z 2 1 Z 3 


A 


-1 


1 


75 


10.5 


16.5 


B 


-1 


1 


40 


8 


14 


C 


-1 


1 


-60 


8 


14 


D 


-1 


1 


-90 


6 


12 
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Figure 14. Near-field source locations. 
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The construction and excitation of the monopole is illustrated in Figure 1 5. Sev- 
eral radiation patterns for the monopole are presented in Appendix A. Computed 
current data is given in Tables 2, 3, 4, and 5 for transmit locations A, B, C, and D 
respectively. 




Figure 15. Transmit monopole and its excitation. 
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TABLE 2. COMPUTED DATA FOR THE ORIGINAL SHIP CONFIGURATION. 

TRANSMIT LOCATION A. 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.81876x1 O' 3 


-71.654 


0.40938 x 1 0' 2 


M-1760 


0.98694 x 10* 4 


-96.254 


0.49347 x 10- 4 


A-106 


0.11971 x IQ' 4 


-39.001 


0.05985 x 1 0' 4 



TABLE 3. COMPUTED DATA FOR THE ORIGINAL SHIP CONFIGURATION. 

TRANSMIT LOCATION B. 



Current Density 
J (A/m ) 


Current 
I (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.19715 x 10- 3 


29.841 


0.098575 x 1 0’ 3 


M-1760 


0.52686 x 1 0' 4 


144.322 


0.26343 x 10- 4 


A-106 


0.58949 x 1 0' 5 


-121.743 


0.294745 x 1 0' 5 



TABLE 4. COMPUTED DATA FOR THE ORIGINAL SHIP CONFIGURATION 



TRANSMIT LOCATION C. 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.22133 x 1 0' 4 


-2.144 


0.1 10665 x 10' 4 


M-1760 


0.72595 x 1 0' 3 


-101.293 


0.362975 x 1 0' 3 


A-106 


0.20377x 1 O' 3 


-59.492 


0.101885 x 10' 3 



TABLE 5. COMPUTED DATA FOR THE ORIGINAL SHIP CONFIGURATION 



TRANSMIT LOCATION D. 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.23681 x 10' 4 


167.264 


0.1 18405 x 10' 4 


M-1760 


0.59780 x 1 0' 3 


63.377 


0.2989 x 1 0' 3 


A-106 


0.93215 x 10' 3 


42.371 


0.466075 x 1 0' 3 
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Next, modification of the ship is simulated by adding a cylinder (possibly 
representing a gun) at either the bow or the stern (see Figures 16 through 18). 
The cylinder was constructed with height of 3m and radius of 1m. The current 
calculations described above are repeated for the same locations. The results 
are summarized in Tables 6, 7, 8, and 9 for the forward modification, and Tables 
10, 1 1, 12, and 13 for the aft modification. 
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Figure 16. Simulated modifications on the topside of the ship ( E and F). 
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Figure 17. Enlarged view of the modification at the bow. 
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Figure 1 8 . Enlarged view of the modification at the stem. 
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TABLE 6. COMPUTED DATA FOR THE MODIFIED (BOW) CONFIGURATION. 



TRANSMIT LOCATION A 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.83021 x 10' 3 


-73.414 


0.41510 x 10- 3 


M-1760 


0.99373 x 1 0’ 4 


-96.343 


0.49686 x 1 0' 4 


A-106 


0.12049 x '\ 0' 4 


-39.060 


0.06024 x 1 0' 4 



TABLE 7. COMPUTED DATA FOR THE MODIFIED (BOW) CONFIGURATION. 

TRANSMIT LOCATION B. 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.19689 x 1 0' 3 


29.080 


0.098445 x 10’ 3 


M-1760 


0.55421 x 10- 4 


146.557 


0.277105 x 10- 4 


A-106 


0.58293 x 10- 5 


-122.617 


0.291465 x 10- 5 



TABLE 8. COMPUTED DATA FOR THE MODIFIED (BOW) CONFIGURATION. 

TRANSMIT LOCATION C. 



Current Density 
J ( A/m) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.22026 x 10' 4 


-2.532 


0.11013 x 10- 4 


M-1760 


0.72598 x 10- 3 


-101.289 


0.36299 x 1 0’ 3 


A-106 


0.20377 x IQ' 3 


-59.489 


0.101885 x 1 0' 3 



TABLE 9. COMPUTED DATA FOR THE MODIFIED (BOW) CONFIGURATION. 

TRANSMIT LOCATION D. 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.23560 x 10- 4 


166.910 


0.1178 x 10- 4 


M-1760 


0.59751 x 10‘ 3 


64.448 


0.298755 x 1 0' 3 


A-106 


0.93214 x 10’ 3 


42.371 


0.46607 x 10‘ 3 
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TABLE 10. COMPUTED DATA FOR THE MODIFIED (AFT) CONFIGURATION. 

TRANSMIT LOCATION A. 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0. 1 0592 x 1 0‘ 2 


-67.518 


0.05296 x 1 0’ 2 


M-1760 


0.13008 x 1 0' 3 


-95.200 


0.06504 x 1 0' 3 


A-106 


0. 1 5443 x 1 0' 4 


-33.833 


0.077215 x 10- 4 



TABLE 1 1 . COMPUTED DATA FOR THE MODIFIED (AFT) CONFIGURATION. 



TRANSMIT LOCATION B. 



Current Density 
J (A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.13313 xIO’ 3 


22.939 


0.066565 x 10' 3 


M-1760 


0.40565x1 O' 4 


144.390 


0.202825 x 10- 4 


A-106 


0.42751 x 10' 5 


-125.496 


0.213755 x 1 0’ 5 



TABLE 12. COMPUTED DATA FOR THE MODIFIED (AFT) CONFIGURATION. 

TRANSMIT LOCATION C. 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.17552 x 10- 4 


-2.643 


0.08776 x 1 0' 4 


M-1760 


0.57529 x 10- 3 


-101.908 


0.287645 x 1 0' 3 


A-106 


0.19817 x 10‘ 3 


-65.225 


0.099085 x 1 0' 3 



TABLE 13. COMPUTED DATA FOR THE MODIFIED (AFT) CONFIGURATION. 

TRANSMIT LOCATION D. 



Current Density 
J (A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.23742 x 10- 4 


167.151 


0.11871 x 10' 4 


M-1760 


0.59929 x 1 0' 3 


64.315 


0.299645 x 1 0' 3 


A-106 


0.92695 x 10’ 3 


42.358 


0.463475 x 10’ 3 
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The difference in the values of the current that were observed can be 



quantized in terms of their fractional change which is defined as follows: 
FRACTIONAL CHANGE = x 100. (3) 

' original 

The calculations for the fractional changes were performed only for the cases of 
the deck-edge antennas at the extremes of the ship, because these are repre- 
sentative of the largest and smallest antenna-source distances. The results are 
summarized in Table 14. 



TABLE 14. FRACTIONAL CHANGES IN THE DECK-EDGE ANTENNA 

CURRENTS. 



DECK-EDGE 

ANTENNA 

LOCATION 


SOURCE 

LOCATION 


MODIFICATION 


FRACTIONAL 

CHANGE 


B-10 


A 


BOW 


1.39% 


B-10 


B 


BOW 


0.13% 


B-10 


C 


BOW 


0.48% 


B-10 


D 


BOW 


0.51% 


A- 106 


A 


ASTERN 


2.9% 


A- 106 


B 


ASTERN 


2.7% 


A- 106 


C 


ASTERN 


2.74% 


A-106 


D 


ASTERN 


0.55% 



The excitation of 100V was chosen to give approximately the same order of 

magnitude power density at the nearest deck-edge antenna. For the current at- 
sea calibration, the transmit site uses a log-periodic dipole (G * 2 dBi) and a 
CW power of P t =500W. The field strength at the ship, with the ship being at a 

distance of R=10nm from the transmitting site (Figure 8), is 
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IEI = ^°g G< = 0.00832 Vim (4) 

and therefore, the power density at the ship is 

P d = = 1 .84x1 0- 7 Wlm 2 . (5) 

For comparison, the power densities at the stern deck-edge antenna (edge 
number A-106) were computed for the near-field excitations A and D. The data is 
shown in Table 15 at the coordinates x=7.2 , y=-84 , z=4.8. Both bow and stern 
excitations with the sources residing at the far ends of the original ship (locations 
A and D) were considered. 



TABLE 15. COMPUTED ELECTRIC AND MAGNETIC FIELD COMPONENTS. 



SOURCE 

LOCATION 


ELECTRIC and MAGNETIC FIELD 
COMPONENTS 


MAGNITUDE of 
ELECTRIC and 
MAGNETIC FIELD 
COMPONENTS 


A 


Ex=0.64648x1 0'-jO. 1 5699x1 O' 2 


|Ex|=0.0647 


Ey=0.521 09x1 0' 2 -j0. 68087x1 O' 3 


|Ey|=0.0053 


Ez=-0. 33703x1 0' 1 +j0.26763x1 O' 2 


|Ez|=0.0038 


Hx=0. 744887x1 O^+jO. 23441 xl O' 4 


|Hx|=7.8x10' 5 


Hy=-0.42164x10' s -j0.20512x10- 4 


|Hy|=2.09x10' 5 


Hz=-0. 1 21 20x1 0' 3 -j0.62355x1 O' 4 


|Hz|=1.36x10' 4 


D 


Ex=-0. 16615-jO. 56306 


|Ex|=0.5871 


Ey=0. 1 5930x1 0' 1 -j0.46265x1 0 3 


|Ey|=0.0159 


Ez=-0. 85361x1 0' 1 -jO. 27030 


|Ez|=0.2835 


Hx=-0.26022x1 0' 3 -j0.47239x1 O' 3 


|Hx|=5.4x10' 4 


Hy=0.71 1 99x1 0' 3 +j0.6741 0x1 O' 3 


|Hy|=9.8x10' 4 


Hz=0. 5598 1 xl 0' 3 +j0. 1 0948x1 O' 2 


|Hz|=0.0012 
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The average power density can be calculated using the formula P av =|E| 2 /r) > 
where r| is the free-space impedance ( r|=377Q ). The power densities due to 
the vertical component of the E-field ( E z ) and the total E-field are shown in 
Table 16. 



TABLE 16. COMPUTED POWER DENSITIES. 



SOURCE LOCATION 


|EJ 2 /ri 


|E| 2 /r, 


A 


3.83023x10-® W/m 2 


3.47x10-® W/m 2 


D 


2. 131 88x1 O' 4 W/m 2 


3.04x1 O' 3 W/m 2 



These densities are roughly comparable to the power densities in the far-field 
test method given by (5), at least at location D. 



C. PLANE WAVE SIMULATION 

The second part of this simulation examines the response of the antenna 
elements due to an incident plane wave. This provides an estimate of the frac- 
tional change in current typical of the at-sea calibration method. Wave angles 
considered range from the bow for 0 ( =90° and <j> =90° , to broadside for 0 t =90° 
and <j>=0°. The measurements were made at the same edges as previously 
noted for both the original and the modified model, and are shown in the follow- 
ing Tables 17 through 25. 
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TABLE 17. COMPUTATION DATA FOR THE ORIGINAL CONFIGURATION. 
PLANE WAVE EXCITATION (0, = 90° , <J>= 0° ). 



Current Density 
J ( A/m) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.82586 x 1 0' 2 


-151.39 


0.41293x1 0' 2 


M-1760 


0.97366 x 1 0' 2 


54.111 


0.48683 x 10' 2 


A- 106 


0.84297 x IQ’ 2 


-122.784 


0.42148 x 1 0' 2 



TABLE 18. COMPUTATION DATA FOR THE ORIGINAL CONFIGURATION. 
PLANE WAVE EXCITATION (0 t = 90° , 4>=45° ). 



Current Density 
J (A/m) 


Current 

1(A) 


Edge 


Magnitude 


Phase 


l=0.5x| J | 


B-10 


0.77574x1 0' 2 


-151.43 


0.38787x1 O' 2 


M-1760 


0.38094x1 O' 2 


-88.6 


0.1 9047x1 O' 2 


A-106 


0.531 85x1 O' 2 


-162.08 


0.26592x1 O' 2 



TABLE 19. COMPUTATION DATA FOR THE ORIGINAL CONFIGURATION. 
PLANE WAVE EXCITATION (0, = 90° , <(j=90 0 ). 



Current Density 
J (A/m) 


Current 

1(A) 


Edge 


Magnitude 


Phase 


l=0.5x|J| 


B-10 


0.45823x1 O' 2 


-106.090 


0.22.91 1x1 0‘ 2 


M-1760 


0.1 8364x1 O' 2 


-43.143 


0.091 82x1 0‘ 2 


A-106 


0.74648x1 O' 3 


-23.478 


0.37324x1 O' 3 
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TABLE 20. COMPUTED DATA FOR THE MODIFIED (BOW) CONFIGURATION. 
PLANE WAVE EXCITATION (0, = 90° , <J> = 0°). 



Current Density 
J (A/m) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.83451 x 10' 2 


-151.03 


0.41725 x 1 0' 2 


M-1760 


0.97341 x 10' 2 


54.117 


0.48670 x 1 0' 2 


A-106 


0.84300 x 10' 2 


-122.783 


0.4215 x 10' 2 



TABLE 21. COMPUTED DATA FOR THE MODIFIED (BOW) CONFIGURATION. 
PLANE WAVE EXCITATION (0 t = 90° , <J>=45° ). 



Current Density Current 

J ( A/m ) 1(A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.7731 4x 10- 2 


-151.34 


0.38657X 10- 2 


M-1760 


0.38088x 10' 2 


-88.59 


0.19044x 10- 2 


A-106 


0.53 185x1 0‘ 2 


-162.08 


0.265925x 10 2 



TABLE 22. COMPUTED DATA FOR THE MODIFIED (BOW) CONFIGURATION. 
PLANE WAVE EXCITATION (0 ( = 90° , <J>=90° ). 





Current 


Density 


Current 




J 


(A/m) 


1 (A/m) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.47561x1 O' 2 


-107.37 


0.23780x1 0' 2 


M-1760 


0.1 8269x1 0' 2 


-43.5 


0.091 34x1 0 2 


A-106 


0.74647x1 O' 3 


-23.62 


0.37323x1 O’ 3 
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TABLE 23. COMPUTED DATA FOR THE MODIFIED (AFT) CONFIGURATION. 
PLANE WAVE EXCITATION (0 i = 90° , <\>= 0° ). 



Current Density Current 

J ( A/m ) 1(A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.82579X 10' 2 


-151.39 


0.41 289x 10- 2 


M-1760 


0.97384X 10' 2 


53.97 


0.48692x 10‘ 2 


A-106 


0.84502X 10' 2 


-122.92 


0.42251 x IQ 2 



TABLE 24. COMPUTED DATA FOR THE MODIFIED (AFT) CONFIGURATION. 
PLANE WAVE EXCITATION (0, = 90° , 4> = 45° ). 



Current Density 
J ( A/m ) 


Current 
1 (A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.77576 x 1 0' 2 


151.423 


0.38788 x 1 0‘ 2 


M-1760 


0.38185 x 1 0' 2 


-88.395 


0. 1 9092 x 1 0' 2 


A-106 


0.53296 x 1 0' 2 


162.172 


0.26648 x 1 0' 2 



TABLE 25. COMPUTED DATA FOR THE MODIFIED (AFT) CONFIGURATION. 
PLANE WAVE EXCITATION (0 ( = 90° , <}> = 90° ). 



Current Density Current 

J ( A/m ) 1(A) 


Edge 


Magnitude 


Phase 


1 = 0.5 x | J | 


B-10 


0.45830 x 1 0' 2 


-106.09 


0.22915 x 10' 2 


M-1760 


0.18182 x 10' 2 


-43.09 


0.09091 x 10‘ 2 


A-106 


0.75617 x 10' 3 


-22.87 


0.37808 x 1 0' 2 
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In order to quantify the severity of the modifications, the fractional changes 
are calculated as in the previous section. Similarly, as in the previous section, 
the calculations were performed only for the deck-edge antennas at the extreme 
front and back of the ship; that is, for the antennas residing on edges B-10 and 
A-106. The results of these calculations are listed in Table 26. 



TABLE 26. FRACTIONAL CHANGES IN THE DECK-EDGE ANTENNA 

CURRENTS. 



DECK-EDGE 


PLANE WAVE 


MODIFICATION 


FRACTIONAL 


ANTENNA 


DIRECTION OF 




CHANGE 


LOCATION 


INCIDENCE 






B-10 


JP 

II 

CD 

O 

o 

je- 

ll 

o 

o 


BOW 


1.05% 


B-10 


0. = 90° <j>. = 0° 


ASTERN 


0.0035% 


B-10 


JP 

II 

CD 

O 

o 

^e- 

n 

cn 

o 


BOW 


0.33% 


B-10 


0. = 90° (j). = 45° 


ASTERN 


0.0025% 


B-10 


9j = 90° 4, = 90° 


BOW 


3.79% 


B-10 


0 i = 90° (J), = 90° 


ASTERN 


1 .29% 


A-106 


0, = 90° h = 0° 


BOW 


0.0035% 


A-106 


0, = 90° 4> j = 0° 


ASTERN 


0.24% 


A-106 


0. = go 0 4>. = 45° 


BOW 


0% 


A-106 


0. = go 0 <J>. = 45° 


ASTERN 


0.20% 


A-106 


_CD 

II 

CD 

O 

o 

je- 

ll 

CD 

9, 


BOW 


0.0013% 


A-106 


0i = 90° <j>, = 90° 


ASTERN 


1 .29% 
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V. CONCLUSIONS 



In this study, an investigation of a new approach in dealing with the "need 
to calibrate" problem of shipboard HFDF antennas was conducted. This ap- 
proach is based on measuring the system response to multiple on-board near- 
field sources. The test is performed along with the standard calibration, and the 
antenna responses to the near-field sources are stored in a reference database. 
Whenever a modification is made to the topside, the near-field test is repeated 
and the new results are compared to the reference database. A significant differ- 
ence may indicate a need to perform a full system calibration. 

In order to validate the proposed approach, a simulation was performed, 
which can be described as a three step procedure: 

1 . compute the response of the antennas due to a source placed at various 
locations on the ship within the near-field of the antennas, 

2. perform the same set of computations after having introduced a simulated 
modification to the topside of the original model, and 

3. compare the results to determine whether any difference is significant. 

The simulation results were characterized by the fractional change in cur- 
rent across selected edges before and after modifications. To provide a refer- 
ence, the at-sea calibration was also simulated and the fractional changes found 
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to be comparable to the near-field case. In fact, the fractional change tended to 
be larger for the near-field excitation than the plane wave excitation. 

An attempt was made to keep the power densities approximately equal for 
the plane wave and near-field excitations. In spite of this, the current densities 
induced on the ship differed by approximately an order of magnitude (higher in 
the plane wave case). 

Based on the results of the few calculations performed in this study, it ap- 
pears that the near-field method is more sensitive to topside changes than the 
at-sea method. This has both advantages and disadvantages. The advantage is 
that if no fractional change is observed for a topside modification, then one can 
be confident that the calibration voltages are still valid. On the other hand, small 
fractional changes in the near-field case can probably be disregarded. However, 
it is crucial that a fractional change "threshold" be established. 

In an effort to set a calibration threshold based on fractional changes in the 
near-field, further comparisons should be made between the near-field and the 
plane wave excitations. Furthermore, both significant and insignificant modifica- 
tions should be examined. A significant modification is one that is known to 
cause the reference voltages to become invalid; an insignificant modification 
does not invalidate the calibration. Other areas that should be investigated are 
convergance of the PATCH model and calculation of the CIDF coefficients for 
near-field sources. 
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APPENDIX 



The radiation patterns of the monopole at various locations are given in this 
Appendix. These patterns were computed in the 4>=90° plane and for the original 
ship configuration only, since the modifications were shown to have no signifi- 
cant impact on the patterns. The patterns for the different source locations are 
shown in Figures A1 to A8. 
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Figure A1 . Original ship. E(0) pattern for the case of excitation at the bow 

(location A). 
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Figure A2. Original ship. E(<j>) pattern for the case of excitation at the bow 

(location A). 
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|E(theta)| in dB 




Figure A3. Original ship. E(0) pattern for the case of excitation at the bow 

(location B). 
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Figure A4. Original ship. E(<j>) pattern for the case of excitation at the bow 

(location B). 
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Figure A5. Original ship. E(0) pattern for the case of excitation at the stern 

(location C). 
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Figure A6. Original ship. E(<J») pattern for the case of excitation at the stern 

(location C). 
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Figure A7. Original ship. E(0) pattern for the case of excitation at the stern 

(location D). 
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Figure A8. Original ship. E(<j>) pattern for the case of excitation at the stern 

(location D). 
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